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The process of stimulated Raman adiabatic passage (STIRAP) provides a possible route for 
the generation of a coherent molecular Bose-Einstein condensate (BEC) from an atomic BEC. We 
analyze this process in a three-dimensional mean-field theory, including atom-atom interactions and 
non-resonant intermediate levels. We find that the process is feasible, but at larger Rabi frequencies 
than anticipated from a crude single-mode lossless analysis, due to two-photon dephasing caused by 
the atomic interactions. We then identify optimal strategies in STIRAP allowing one to maintain 
high conversion efficiencies with smaller Rabi frequencies and under experimentally less demanding 
conditions. 
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I. INTRODUCTION 

Coherent conversion of an atomic to a molecular Bose- 
Einstein condensate (BEC) is a first step towards 'su- 
perchemistry' ||l||, which is the stimulated emission of 
molecules in a chemical reaction. A number of studies 
of this H, H have shown that direct conversion via 
Raman photoassociation |||] appears feasible, based on 
stimulated free-bound and bound-bound transitions in 
the presence of two laser fields of different frequencies 
Here pairs of atoms from the two-atom continuum of the 
ground electronic potential are transferred - via an ex- 
cited bound molecular state - to a bound molecular state 
of a lower energy in the ground potential. Raman pho- 
toassociation allows coupHng to a single molecular state, 
which can be selected by the Raman laser frequencies. 
Practical estimates using available lasers and transitions 
indicate that coherent transfer may be limited by sponta- 
neous emission from the intermediate molecular excited 
electronic state. Another mechanism that can result in 
coupled atomic-molecular BEC systems |0| is based on 
Feshbach resonances 1^. However, realistic analysis and 
experimental implementations |^ indicate that the loss 
processes due to inelastic atom-molecule collisions occur 
at a significant rate. 

A possible route towards minimizing losses and deco- 
herence from spontaneous emission in photoassociation is 
stimulated Raman adiabatic passage (STIRAP) in 
which a counter-intuitive pulse sequence is used, where 
the first input pulse couples the molecular levels - even 
when there are no molecules present. In this situation, a 
dark superposition state is formed, due to interference ef- 
fects between the atomic and molecular electronic ground 
states. This minimizes the probability of a real transi- 
tion to the molecular excited state, and hence reduces 
spontaneous emission. Previous analyses of this problem 
have not taken into account losses, collisions, or the full 



three-dimensional structure of the two Bose condensates 
in a trap. 

In this paper, we provide an analysis which is rele- 
vant in a physically appropriate model that does include 
the known physics of spontaneous emission losses, s-wave 
scattering processes and spatial diffusion of the conden- 
sates. The result is that the STIRAP process appears 
feasible at high laser pulse intensities, provided the Rabi 
frequency is much greater than the two-photon detuning 
due to mean-field interactions. We give a detailed numer- 
ical calculation based on both homogeneous and inhomo- 
geneous three-dimensional mean-field (Gross-Pitaevskii 
type) theories, including couplings to non-resonant inter- 
mediate levels, and show how the results scale with the 
two-photon detuning, pulse duration and Raman pulse 
intensities. An optimal situation is found by considering 
an off-resonance operation and different effective Rabi 
frequencies in the two Raman channels. We show that 
these strategies can greatly enhance the conversion effi- 
ciency for given laser intensities, thus making the exper- 
imental requirements much more feasible. 



II. COUPLED GROSS-PITAEVSKII 
EQUATIONS FOR STIRAP 

We start by considering the theory of coherently in- 
teracting atomic and molecular condensates needed to 
describe this process [Q, ^, and assume a specific cou- 
pling mechanism based on stimulated free-bound Raman 
transitions |^, in which two atoms of energy Ei collide 
to from a molecule of energy E2 , with an excited molec- 
ular state forming as an intermediate step. The Raman 
coupling is induced by two laser fields of frequencies ui 
and u)2, and becomes resonant when the Raman detun- 
ing S — {2Ei — E2)/h~ {ijj2 — t^i) goes to zero. This 
allows coupling to a single molecular state, which can be 
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terms in the effective interaction Hamiltonian: 
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Figure 1: Diagrammatic representation of the free-bound and 
bound-bound transitions in STIRAP. 



selected by the Raman laser frequencies. 

We derive the atom-molecule coupling for a simplified 
model of the two-body interaction 0, |ll|, in which the 
atoms interact in their electronic ground state through 
a single Born-Oppenheimer potential Vg{R). Molecules 
are formed in a single bound vibrational state of energy 
E2 with radial wave function U2{R). Two free atoms 
with zero relative kinetic energy have a total energy 2Ei , 
and a relative radial wave function ui{R), normalized 
so that asymptotically ui oc {1 — ai/R). We assume 
that the laser field has two frequency components, with 
E = ReX; [e(') exp(iw,t)] , i = 1,2. Each couples the 
ground electronic state to a single electronically excited 
state described by a potential Ve{R), with 'bare' elec- 
tronic Rabi frequencies — d3i(R)-E(')/^ ■• where 
d3,i(R) is the molecular electric dipole matrix element 
with a nuclear separation of R. The excited state has 
vibrational levels \v') with energies E^i and radial wave 
functions Uv' (R) . All bound levels are normalized so that 

We consider the case of near resonant transitions 1 — > 
— > 2 and denote the resonant excited vibrational level 
v' via index 3 (see Fig. [l]). The usual quantum field 
theory Hamiltonian ||T^ for noninteracting atomic (i — 1) 
and molecular {i = 2, 3) species, in well-defined internal 
states described by annihilation operators ^i, is given 
by: 
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|V#,(x)|2-^y,(x)#l(x)^',(x) 



(1) 

Here (i — 1,2,3) are the masses, m2,3 — 2mi, and 
T^i(x) is the trapping potential including the internal en- 
ergy for the i-ih species, where we define Vi{Q) = Ei. 

Including s-wave scattering processes and laser in- 
duced particle inter-conversion, together with the as- 
sumption of a momentum cut-off, results in the following 
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e — *f(x)«'T(x) +i?.c. 
e~^'^^**2(x)'I'5(x) + i/.c. 



(3) 
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Here = ^ (Plinf\ll)ul{R)ui{R) « f^f'^/j,3 {i = 
1,2) are the molecular Rabi frequencies. These can 
be treated using the Franck-Condon overlap integrals 

= J d^'Rul{R)ui{R) , if we take Sl-*^'^ as the mean 
electronic Rabi frequency. We note that fli, which con- 
nects the atomic and molecular condensates, has units 
of s~^TO~^/^, and must be multipHed by the atomic 
condensate amplitude to obtain a true Rabi frequency. 
The couplings Uij in the diagonal case are given by 
Uii — ATThai/rrii, where a.i is the i-th species s-wave scat- 
tering length, while the non-diagonal terms are given by 



Uji = 2'Khaij / fj,ij , where aij is the inter-species 



scattering length and fi = mirrij / [mi -\- rrij) is the re- 
duced mass. 

In addition, we account for losses from each state, at a 
rate 7^. The resulting Heisenberg equations for the field 
operators are treated within the mean-field approxima- 
tion, in which the operators are replaced by their mean 
values, and a factorization is assumed. This approxima- 
tion is expected to be valid at sufficiently high density. 
Corrections due to quantum correlations have been 
treated in greater detail elsewhere |Q. Next, we intro- 
duce rotating frame detunings, defined so that: 



2Ai(x) = {E^~2Vi{^))lh~uoi 
A2(x) = (S3-y2(x))A-c^2 , 
A3(x) - (^3 -V^3(x))A . 



(5) 



In the case of uniform condenstates, Vi^x) are equal to 
V,{x)=E„&nA A3 = 0. 

This results in the following set of Gross-Pitaevskii 
type of equations for the mean-field amplitudes, in ro- 
tating frames such that -01 = ( ) exp[i(£'i -I- Ai(0))i/ft]: 



(6) 



dt 



lA^^02 + ^03, 



= ^A3 V^3 + ^0, + — ^2, 

Here A^^ is the i -th Gross-Pitaevskii mean-field detun- 
ing in the rotating frame, defined so that: 



Af^(x,i) = A,(x)-f ^V2-^t/,fc|^,p-fz| . (7) 
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We also introduce the two-photon laser detuning at trap 
center: 



d = A2(0) - 2Ai(0) = -(£^2 - 2Ei)/h+iLUi - 1^2) ■ (8) 



In addition to losses due to spontaneous emission from 
the electronic excited states, rotationally or vibrationally 
inelastic atom-molecule collisions may also give rise to 
losses. The magnitude of these rates is presently un- 
known, and we neglect them here. We note that these 
rates should decrease rapidly with increasing molecular 
binding energy and go to zero in the molecular ground 
state, so that it should be possible to obtain a very low 
rate by selecting the coupHng to a deeply bound molec- 
ular level. This approximation means that we will set 
7i = jSsi in the following treatment, where 6ij is the 
Kronecker delta-function. 

The simplest STIRAP scheme employs exact tuning of 
the laser frequencies to the bare state resonances, i.e., 
Ai(0) — A2(0) — 0. This, however, is not necessarily 
required as the STIRAP can still occur with detuned in- 
termediate levels. Moreover, as we show below, in BEC 
environments, the off-resonance operation turns out to be 
more efficient if the detunings compensate for the phase 
shifts due to mean-field energies. In effect this is equiv- 
alent to a renormalized two-photon "on-resonant" opera- 
tion in which A^^(x, 0) - 2Af ^(x, 0) « 0. 

If, instead of considering near-resonance coupHng, we 
consider a large intermediate level detuning so that the 
excited state can be adiabatically eliminated, we recover 
the basic terms in the set of equations analyzed in |Q. 
Simultaneously, this would give us the previously known 
result jl^ for the laser induced modification to the scat- 
tering length ai occurring in a simple single-laser pho- 
toassociation of pairs of atoms, thus justifying the above 
form of the interaction Hamiltonian. In the present pa- 
per, however, we assume that the detunings Ai are small 
compared to the characteristic separation between the 
vibrational levels so that all other vibrational levels can 
be neglected. The near-resonant excited level is treated 
explicitly, rather than eliminated adiabatically as in [Q. 



tions: 



III. STIRAP IN A BEC 



Before carrying out simulations of the full 3-D equa- 
tions in a trap, it is instructive to start with a simplified 
version of the theory - in which there are no kinetic en- 
ergy terms. We expect this approximation to be valid 
in the Thomas-Fermi limit of large, relatively dense con- 
densates, which is a regime of much current experimental 
interest. This is described by the following set of equa- 
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dtp3{x,t) 
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= zA^^^2 + ^^3, 



(9) 



where we have introduced as an effective Thomas-Fermi 
limit frequency shift. 



Aj^(x,t) = A,(x)-^C/,,|7^fe 



'+^^- (10) 
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This corresponds to the treatment given in [g^, except 
that we explicitly include the loss term 7 due to spon- 
taneous emission, the s-wave scattering processes due to 
Uij, and the Franck-Condon integrals into the coupling 
coefficients for the free-bound and bound-bound transi- 
tions. 

In order to understand how this is related to the usual 
STIRAP technique in a three-level A atomic system, we 
introduce a new wave-function ip„i = "01 / V2, which cor- 
responds to the coherent amplitude of a (virtual) molec- 
ular condensate with the same number of atoms as in the 
atomic BEC. We then introduce a Bose-stimulated Rabi 
frequency, which includes a local coherent BEC ampli- 
tude for the first free-bound transition: fii = ipl^i. This 
leads to the equations: 
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These are precisely the usual STIRAP equations, except 
with additional detunings coming from the mean-field in- 
teractions, and a Rabi frequency in the first transition 
which is proportional to the amplitude of the atomic BEC 
wave-function. In practise, the Rabi frequency may have 
an additional space-dependence due to the spatial varia- 
tion of the laser phase and intensity. We therefore con- 
clude that, provided we can satisfy the normal adiabatic 
STIRAP requirements of slow time-variation in the ef- 
fective Rabi frequencies in the above equations, the tech- 
nique will also work for a BEC. This is a simpler proof 
than previously |0. In particular, we can immediately 
deduce the expected solution for real Rabi frequencies in 
the adiabatic Hmit. In order to have ^iipm + f^2'02 = 0, 
one requires: 



tpi{x,t) = '0i(x,O)cos(6'), 
M^,t) = -0i(x,O)sin(0)/V2, 
^3(x,i) = 0. 



(12) 
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Here the space dependent mixing angle 0(x, t) is ob- 
tained from the ratio of effective Rabi frequencies: 



tan(^) . - 



no 



n 1/2 



We can see that ini 



(13) 

while ipi is still close to its 



initial value, the mixing angle is close to its expected 
value in normal STIRAP, since tan(6') « 0)ili/r22. 
However, at the final stages of the adiabatic passage, the 
nonlinear effects due to the atom-molecular coupling be- 
come important. As the atomic EEC amplitude only 
varies on the time-scale of the input fields, the nonlinear 
atom-molecular coupling term by itself should not intro- 
duce new adiabatic restrictions. However, there will be 
time-dependent detunings introduced by the mean-field 
terms. 

Finally, we can see that similar conclusions can also be 
reached in a non-uniform BEC, by replacing the uniform 
detunings with appropriate Gross-Pitaevskii detunings, 
that include the spatial potentials. In a typical BEC 
cooled in the Thomas-Fermi regime, we expect the kinetic 
energy terms to have relatively small effects. 

STIRAP can therefore be implemented as usually by 
using two laser pulses applied in counterintuitive order. 
We choose Gaussian pulses of the form 



nf\t) = nr'"'^MHt-ur/T% (z = i,2), (14) 



or 



exp[ 



(15) 



(0) 



where the peak values are related as follows: 

r2^^''°^/i_3. The pulse at frequency uj2 is applied first, with 
the center at t2, while the second pulse at frequency uJi 
is delayed by aT, i. e.. 



ti — aT, 



(16) 



where a is the delay coefficient, and T is the pulse dura- 
tion, which we assume is the same for both pulses. 

In terms of the Rabi frequencies, the adiabatic condi- 
tion for STIRAP now reads as [lisl: 



17(i)AT> J I 



AriAr, 



(17) 



where n{t) = ^|r2i(i)|2 |172(0P is the rms Rabi fre- 
quency. At is the duration during which the pulses over- 
lap, and {2Aj^ — A|^-^) simply corresponds to the de- 
tuning of the single-photon transition. However, there is 
a second condition, which is often not stated expHcitly. 
This is that STIRAP requires an effective two-photon res- 
onance, to avoid dephasing between the initial and final 
states in the dark-state superposition. The two-photon 
resonance condition is different from the usual STIRAP 
case, since a detuning of A^^ causes a phase rotation 




Figure 2: The Rabi frequencies fll^^^^t) for the Raman tran 



sitions, with the peak values of n[''^^'°^ = = 2.1 x 10^ 

s~^, pulse durations T = 10~^ s, and a delay coeflEicient of 
a = 1.5. 



both in 1pm and in fli as well, since this also includes a 
phase term from the condensate. As a result, the nec- 
essary condition for two-photon resonance is therefore: 



.TF 



2Af^|Ar « 1 



(18) 



This leads to a third condition, which shows that there 
is a lower bound to the allowed Rabi frequency in order 
to have STIRAP occurring in the presence of mean-field 
dephasing effects: 



n{t) » |A^^ 



2A^ 



TF\ 



(19) 



As is usually the case in STIRAP, these conditions 
cannot be satisfied very early or late in the pulse se- 
quence, when the Rabi frequencies are small; but they 
should be satisfied over most of the STIRAP interaction, 
and over most of the condensate volume. For simplic- 
ity, we will apply these conditions to the peak Rabi fre- 
quency and to the total pulse duration T. Further, 
since ili is itself a function of the STIRAP evolution, 
we introduce an effective first Rabi frequency, defined in 
terms of the initial density ni(0) = 1-01(0)1^. This is 
an upper bound to the stimulated Rabi frequency; thus 
'^^^^\t) = ^ni(0)17i(t) > hi{t) (and sometimes we 
^'^ffh-f\ — fl2{t), for uniformity), where 



n 



write n'^^'^'it) 



^'"\t) = nl'"-"^ exp[-(t - U)yT% (^ = 1, 2) 



(e//,0) 



f^i'''°^/i.3y^ and 



(20) 



where n\ 

"2 ^2,3- 

A typical pulse sequence is shown in Fig. ^. From the 
definition of fi^, the Franck-Condon overlap integrals are 
an important issue from the point of view of employing 
a realistic set of parameters. Since the overlap integrals 
enter into the definition of the effective Rabi frequencies, 
their values will affect the adiabatic condition rewritten 
in terms of the "bare" electronic Rabi frequencies fi^ . 
We will analyze this in more detail in the next section. 
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7.4 X lO'^ s"^ 


ni(0) 


4.3 X 10^° m"^ 




2.1 X lO'^ s"^ 


T 


10-^ s 



Table I: Typical parameter values for efScient STIRAP. 
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4.96 X 10"" mVs 




-6.44 X 10"" mVs 


U22 


2.48 X 10"" mVs 


(73, 






Table II: Typical mean-field interaction potentials in Rubid- 
ium condensates. 



IV. UNIFORM CONDENSATE RESULTS 

We start by considering a uniform condensate, de- 
scribed by a similar equation to the Thomas-Fermi case, 
except with a uniform trap potential for simplicity: 



dt 

^V'2(x,t) 

dt 

^-03(x,t) 

dt 



in 



(21) 
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Here the uniform detuning term A^^(O) is defined as 
the Thomas-Fermi detuning, evaluated at the trap cen- 
ter. We start by considering a uniform condensate in 
which the s-wave scattering interactions are negligible 
(i.e., Uij = 0), and assume exact resonances with respect 
to bare state transitions, Ai = A2 = 0. We first sim- 
ulate the above simplified equations ( ^ ) with an initial 
condition of a pure atomic condensate. This model is 
more realistic than that of M. Mackie et al. [0, as it 
includes spontaneous emission. We find that including 
the loss term 7 imposes restrictions on the effective Rabi 
frequencies and the duration of the pulses T. 

The results are best understood with reference to Table 
|[ which gives the values of typical STIRAP parameters 
characteristic of a condensate of *^Rb atoms J^, |^), and 
corresponding to the pulse sequence in Fig. gT^ 

Taking the values of the parameters in Table |, and an 
optimum delay coefficient of a ~ 1.5, gives 77 ~ 0.96 or 
about 96% efficiency of conversion of atomic BEC into 
molecular BEC, even including the upper level sponta- 
neous emission. Here, the conversion efficiency 77 is de- 
fined as the fraction of the initial number of atoms ni{0) 
converted into molecules 



V 



2^2(00) 
ni(0) 



(22) 



where ^2(00) is the final number of molecules. This 
accounts for the fact that ni atoms can produce ni/2 
molecules at best. 

For comparison, using smaller Rabi frequencies, 
^(e//,o) ^ 2.1 X 10^ s-\ and a larger value of T = IQ-^ s, 
so that the product fl['^^'^'^''T still has the previous value 
^(eff,o)rp _ 2.1 X 10^, gives a maximum conversion effi- 
ciency of ?7 ~ 0.83, with a new optimum delay coefficient 
a ~ 1.2. This is smaller than in the previous example. 



In order to reach the same efficiency as before, one has 
to further increase the pulse durations (up to T = 10"^ 
s), i.e. enter into a deeper adiabatic regime. In in the ab- 
sence of the spontaneous emmission term, the conversion 
efficiencies would not be different in these two examples. 

In other words, in this simplified model it is possible 
to have effective Rabi frequencies smaller than the spon- 
taneous emission rate 7, provided the duration of the 
pulses is long enough. As usual, we can understand this 
physically as implying that the upper level is never ac- 
tually occupied for very slowly varying adiabatic pulses. 
Hence, just as in the case of atomic STIRAP, we can 
ignore spontaneous emission from the upper level pro- 
vided that we use very slowly varying pulses which are 
sufficiently deep in the adiabatic Hmit. As we see in the 
following calculations, the problem with this strategy is 
that very long pulses will tend to cause violations of the 
two-photon resonance condition, in the presence of mean- 
field interactions. 



A. Effects of the mean-field energies 

We now wish to include the mean-field energy terms, 
and first restrict our analysis to the atom-atom scattering 
processes. We consider a characteristic value of Un — 
4.96 X 10"^^ m'Vs corresponding to the scattering length 
of ^^Rb atoms ||6| oi = 5.4 nm (mi = 1.443 x 10"^^ kg). 
Together with the choice of the initial atomic density 
ni(0) as before (see Table |), the value of Un sets up a 
characteristic dephasing time scale 



tph = [Uiini{0)Y 



(23) 



equal in this case to tph — 4.7 x 10^^ s. The pulse du- 
ration T must be smaller than or of the order of the 
dephasing time, in order to permit STIRAP, otherwise 
the two-photon resonance condition will not be satisfied. 

Thus, including atom- atom scattering imposes an up- 
per limit to the pulse duration, so that T has now to 
satisfy T < tph — 10^^ — lO^"' s. But this restriction 
means that one can no longer use longer pulse durations 
for small values of while still maintaining high 

conversion efficiency. As a result, the adiabatic condition 



(e//,0) 



T > 1, with the restriction of T < 10"^ - 10" 



s, requires high peak values of the effective Rabi frequen- 
cies: > 10^ s-i. 

In order to satisfy this combination of requirements, 
we use typical parameter values given in Table |, with 
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Figure 3: Efficient conversion of an atomic condensate into a 
molecular condensate during STIRAP as obtained by sim- 
ulating Eqs. (p l|), with an initial atom number density 
ni(0) = 4.3 X 10^^ m"^. Other parameter values are as in 
Fig. 2, and Tables § and ^ The solid line indicates atomic 
density, the dashed line the molecular density, and the dotted 
line the analytic result in the adiabatic limit. 



Uii = 4.96 X lO^i"^ m^/s. Simulating Eqs. (|2|) with 
these parameter values and with all other couplings Uij 
set to zero, gives a maximum of ~ 0.95 conversion 
efRciency, for the optimum delay coefficient of a ~ 1.5. 

As the next step, one can include the mean-field en- 
ergies due to atom-molecule {U12, and molecule- 
molecule {U22, U23, U33) scattering processes. Consider- 
ing the fact that the excited molecular state never gets 
highly populated in STIRAP, only processes described by 
U12 and U22 are to be taken into account here. Provided 
that the scattering lengths for these processes are of the 
same order of magnitude as the atom-atom scattering 
length, these terms do not lead to a dramatic change in 
the conversion efficiency. 

To account for the most recent experimental data on 
ultracold atom-molecule scattering in a ®''Rb condensate 
1^, we have included the U12 term with 012 = —9.346 nm. 
In addition, we include the U22 term with an assumption 
that 02 = ai and set U^i = since these are currently 
not known. The resulting values of Uij are summarized 
in Table 0. 

The results of simulations are given in Fig. ||, where we 
see about 93% (77 ~ 0.93) conversion of the atomic con- 
densate into the condensate of molecules, for fl^-^^-^'^^ — 
^(e//,o) ^ 2.1 X 10^ s-i, T = 10-4 s, and a ~ 1.5. This 
figure also includes the analytic theory calculated in the 
adiabatic limit for comparison, and shows that for these 
parameters, the results of the numerical simulation are 
close to those from the adiabatic theory. 

Thus, we conclude that even including the mean field 
energies STIRAP can be carried out, provided one uses 
faster time scales than in the absence of the s-wave scat- 
tering. As a consequence the effective Rabi frequencies 
have to be kept at a rather high value. Characteristic re- 
sults for comparison are summarized in Fig. |^, where we 
plot conversion efRciency ?y versus the relative delay co- 
efficient a, for cases where s-wave scattering are present 




Figure 4: The conversion eflEciency 77 as a function of relative 
delay a, for: (a) T = 10"* s and Uij = 0; (b) T = 10"^ s 
and Uij = 0; (c) T = 10"" s and Uj as in Table @; (d) 
T = 10"^ s and Uj as in Table 0. The full, dashed, and 
dashed-dotted lines correspond to effective Rabi frequencies 
nj"^^'"^ = n'-^^^'"^ equal to 2.1 x 10'' s-\ 6.3 x lO" s-\ and 
2.1 X 10^ s"^. Other parameter values are as in Table |. 



or absent, and for different values of the effective Rabi 
frequencies and pulse durations T. 

Figure || (a) shows a reasonably efficient conversion, 
in the absence of mean-field interactions, but including 
losses. As expected, spontaneous emission losses are re- 
duced, and efRciency is improved further by the use of 
longer pulses, further into the adiabatic limit, as in Fig. 
4 (b). However, the more realistic example given in Fig. 
i (c) which includes mean-field interactions shows rather 
poor conversion, especially when we use smaller effective 



Rabi frequencies, 



(e//,0) 



(e//,0) 



= 2.1 X lO*' 



shown by the full line where maximum ?; ~ 0.12 at opti- 
mum a ~ 0.7. This is caused by the effective two-photon 
detunings induced by the mean-field interactions. The 
situation is made even worse rather than better in Fig. ^ 
(d) when longer pulses are chosen, giving more time for 
two-photon detunings to occur. 

To be more specific about values of the effective Rabi 
frequencies we recall that the definition of fl^^^^^ involves 
the Franck-Condon overlap integrals /^^a and "bare" elec- 
tronic Rabi frequencies = \dM ■ E,;|/?i. Given the 
values of cIm and /^^s which are specific for particu- 
lar dimer species involved, the size of fl'f'^'^'^'^ can be 
translated to the intensities of the Raman lasers. Con- 
sidering ^^Rb2 as an example, and using a characteris- 
tic values of the corresponding Franck-Condon integrals, 
|/i,3| ~ 10-1"* m3/2 and - 0.1 0, the magnitudes of 



^(e//,0) 



= 2.1 X 10 s ^ translate to peak val- 



ues of the "bare" Rabi frequencies equal to Vl^^^''^^ = 10^^ 
s-i (for ni(0) = 4.3 x 10^° m'^) and d^^'^'^ = 2.1 x 10« 



The peak Rabi frequency of ri^*^''"-* 



lOi 



for 



7 



0.4 



0.3 




5 10 

-5 (10^^ s"^) 

Figure 5: The conversion efficiency as a function of —5, 
for q['"'°'> = ni'"'°^ = 2.1 X 10« s-\ T = 10-^ s and 
different delay coefficients: a = 0.2 (full line), a = 0.4 (dashed 
line), a = 0.8 (dashed-dotted line), a — 1.5 (dots). Other 
parameter values are as in Tables H and E3. 



the free-bound transition would be realized with a 1 W 
laser power and a waist size of about 10 /im, which is 
not impossible - but much higher than we would esti- 
mate without the combined effects of spontaneous emis- 
sion and collisional processes. Another obvious problem 
here is that the waist size of 10 /im is comparable to 
characteristic spatial extend of current BECs in a trap. 

In summary, our analysis shows that the relatively 
small overlap integrals for the free-bound transitions, to- 
gether with the mean-field interaction detunings, can re- 
quire rather high intensity of the wi-laser for obtaining 
high conversion efficiencies. 



ergy is changed dynamically as the populations of the 
atomic and the molecular condensate themselves are be- 
ing changed during STIRAP. As a crude estimate of an 
appropriate value of S one can simply choose it to com- 
pensate the initial mean field energy in the atomic con- 
densate. This approach - employed for smaller Rabi 
frequencies than before - can substantially improve the 
conversion efficiency, compared to the case of zero two- 
photon detuning. 

To show this we have carried out simulations with ten 
times smaller Rabi frequencies than before (i.e. with 
^(e//,o) _ 2 1 X 10^ s ^^), corresponding to a decrease 
of the Raman laser intensities by a factor of 100. The 
results are summarized in Fig. || , where we plot the con- 
version efficiency rj versus S, for T — 10^^ s and different 
delay coefficients a. 

As we see, by varying the two-photon detuning and 
tuning it to the optimum value one can improve the con- 
version efficiency by about a factor of two or more, for a 
range of values of the delay coefficient. Furthermore, as 
the effective dephasing time is increased when the con- 
tribution of the mean field energies is compensated by S, 
one can further improve the results by employing longer 
pulse durations. More generally, the problem of finding a 
set of values of T, a, and S that maximize the conversion 
efficiency, for given values of the effective Rabi frequen- 
cies, is now transformed to an optimization problem. 



C. Asymmetric effective Rabi frequencies 



B. Off-resonance operation 

In order to allow one to operate under less demanding 
laser powers or smaller Rabi frequencies (e.g. il'f^^'^^ — 
^(^ff'O) = 2.1 X 10^ s~^) - while still maintaining efficient 
conversion - we now consider the role of the detunings Ai 
and A2 in the off-resonance regime of operation. In effect, 
this approach relies on compensating for the phase shifts 
due to the mean field energies, and tuning the free-bound 
and bound-bound transitions to a "true" resonance. The 
physics behind this is that in BEC environments it is not 
appropriate to consider transitions with respect to single- 
particle bare energies Ei. Rather, the relevant energies 
and therefore the effective resonances have to take into 
account the mean-field energy contributions due to self- 
and cross- interactions between the condensates. 

More specifically, it is the two-photon detuning 5 that 
has to be adjusted to the relative phase between the 
atomic and molecular condensates. Alternatively speak- 
ing, by tuning the two-photon detuning to compensate 
for the net mean-field energy, one reduces the effect of 
dephasing since the effective dephasing time becomes 
longer compared to the pulse durations. The problem, 
however, is more complicated because the mean field en- 



We now wish to explore an alternative strategy for im- 
proving the conversion efficiency under experimentally 
less demanding conditions of smaller Rabi frequencies. 
We consider the effects of non-equal effective Rabi fre- 
quencies. 

Using the earlier given characteristic values of the 
Franck-Condon overlap integrals, |/i,3| ~ lO'^^m^/^ and 
1/2,3! — 0-1 ) ■^fi can estimate that the moderate magni- 
tudes of n[''^^'°'' = f}^"^-^'"' = 2.1 X 106 s-i translate 
to the following peak values of the bare Rabi frequen- 
cies: ni''^'^^ = IQio s-i (for ni(0) - 4.3 x lO^o m'^) and 
^(e;,o) _ 2 2 X 10'' s^^. As we see, while this corresponds 
to equal effective Rabi frequencies, however, the absolute 
values of the corresponding bare Rabi frequencies are not 
equal. The limitation on laser intensities refers primarily 
to the free-bound transition, whose bare Rabi frequency 

As far as the second Rabi frequency fl''^^'^'' is con- 
cerned, one can in principle increase its magnitude up 
to the same value as i.e. ^1^2'^'^'^ = 10^° s~^, thus 

maintaining experimentally similar and reasonably high 
intensities for both lasers. Under these conditions, and 
for the same values of the Franck-Condon overlap inte- 
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Figure 6: The maximum conversion efficiency rj versus Q2 
(the evaluated points are represented by circles), for Q^^'''^^ = 
10^° s~^ and the cor resp onding optimum values of T, a, and 
as given m Table Bill The triangle gives the result of an 
optimization with S — 0. 



is-') 


r (xio-^ s) 


a 


5 (xlO* s"') 


V 


1.5 X lO'^ 


0.987 


0.753 


3.58 


0.257 


2.1 X lO'^ 


0.966 


0.795 


3.57 


0.331 


3 X lO'^ 


1.05 


0.882 


3.44 


0.391 


5 X lO'^ 


1.29 


1.05 


3.11 


0.437 


7.5 X 10'^ 


1.49 


1.19 


2.92 


0.453 


10» 


1.61 


1.30 


2.83 


0.459 


2 X 10* 


1.86 


1.53 


2.71 


0.467 


10^ 


2.32 


1.98 


2.59 


0.474 


10"' 


2.58 


2.51 


2.93 


0.486 



Table III: Optimum STIRAP parameters for: = 
10^°s-\ I/1.3I = 10-" m^^Vs, and ni(0) = 4.3 x lO^" m-^ 
so that rij^'^^^'"' = 2.1 X lO^s"^ in all cases; different values 
of n'^"'^^ = A.sfi^"''"' are taken for |/i,3| = 0.1 and ni"''"'' 
ranging from 1.5 x 10^ to 10^'' s~^. 



grals and ni(0), we would have 



REALISTIC CONDENSATE MODELS 



(ei,0) 



{eff.O) 



2.1 X 10*^s^\ 



(24) 



We can now ask the question of what happens in 
STIRAP with different effective Rabi frequencies, and 
whether one can achieve higher conversion efficiencies in 
the regime where 0,^2^-^^'°^ > This approach 

again leads to an increased conversion efficiency com- 
pared to the case of equal effective Rabi frequencies. To 
generalize the analysis, we now treat different cases as 
an optimization problem (that maximize ry), carried out 
with respect to T, a, and S, for a set of different values 



of O 



(el,0) 



within a range of fl 



(el,0) 



= 1.5 X 10^ - IQi 



110 



In 



s ^, and for a given value of ' — 10 
terms of the effective Rabi frequencies, this corresponds 
to Vl^^^^ ranging from 1.5 x 10^ to 10^ s"^, for a given 

^(e//,0) ^ 2.1 X 106 S-1. 

The results are summarized in Fig. ^ and in Table 01 
where we give the corresponding optimum values of T, a, 
and (5, and the resulting maximum conversion efficiency 

For comparison, in the symmetric case of ^'f^^''^^ = 
^i^fffi) ^ 2.1 X 10^ s-i and (5 = 0, such that the op- 
timization is carried out only with respect to T and a, 
the maximum conversion efficiency would be 77 ~ 0.14 (at 
optimum T = 0.46 x lO""' s and a = 0.54). This case is 
represented by the triangle, in Fig. ||. 

Thus, we have shown that by introducing the pos- 
sibility of varying the two-photon detuning 8 and the 
Rabi frequency fij'^''^'', the conversion efficiency can be 
increased almost by a factor of 4. This can be crucial for 
experimental observation of the phenomenon of coherent 
conversion of an atomic BEC into a molecular BEC, via 
STIRAP. 



A. Uniform multi-level model 

In our model for STIRAP we only treated the coupHng 
of laser lui to the free-bound transition |1) ^ |3) with 
Rabi frequency Vti = ^\ /i 3 together with the cou- 
pling of laser W2 to the bound-bound transition |2) ^ |3) 

- (el) 

with Rabi frequency ^2 = ^2 -^2,3- This approximation 
can only be valid if the laser uji is far detuned from the 
1 2) 4-> 1 3) transition, and similarly - if the laser u)2 is far 
detuned from the |1) ^ |3) transition. In addition, the 
two lasers have to be far detuned from transitions to any 
other vibrational levels \v') (adjacent to |3)) in the ex- 
cited potential. We define the relevant detunings, for the 
simplest uniform follows: 



Ai3,c.. = {E:i-2Ei)/h~L02, 

A23,c.i = [Es- E2)/h-u;i, 

Ai«',u.i = (E^' -2Ei)/h-u;i, 

^iv',u>._ = (E^' -2Ei)/h-u;2, 

^2v',uJi = {Ev' - E2)/fl- UJi, 

^2v',i02 = {E^' - E2)/h- UJ2- 



(25) 



In general, these cross-couplings - if included into the 
model - lead to incoherent radiative losses of atoms and 
molecules due to spontaneous emission, which modifies 
the effective detunings to: 



A2 



ij/2. 



(26) 



In order that these losses be neghgible we require the re- 
spective detunings to be large enough. This requirement, 
however, may not be easily satisfied, as the magnitudes of 
the detunings are in principle limited from above by the 
characteristic distance between the adjacent vibrational 
levels of the excited molecular potential. In other words, 
increasing the detuning with respect to one transition 
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will eventually bring the laser frequency to a resonance 
with respect to the nearby level. More importantly, these 
cross-coupHngs provide scattering pathways that are not 
cancelled out in a dark-state interference effect, so that 
their overall disruptive effect - over the adiabatically long 
pulse durations - may turn out to be rather large. 

In order to estimate these effects, we therefore explic- 
itly include all other relevant coupling processes into our 
model. In addition to losses, the incoherent couplings 
induce Hght shifts that effectively lead to a dephasing 
between the atomic and molecular condensates. Treat- 
ing these, leads to the following additional terms in the 
STIRAP equations, in same rotating frames as in Eqs. 



dt 



■■^"^''^^'''^%-''^"Vi>3, (27) 



— = {...)+ t(3]ij2 



V2 



.,|^2 + ,l^i;^e--*^3, (28) 



dt - 2V2 ' 



(29) 



where (...) stand for the terms already present in the 
right hand sides of Eqs. (^, and oju = lui — uj2, ■ Here, 
the effective complex light shift coefficients (3^ and /?2 , 
the nonlinear phase shift Uu (which effectively leads to a 
modified atom-atom scattering length) and the effective 
parametric coupHngs x, %' (including only non-oscillating 
terms) are given by: 



/37 



(A) 12 



(^)|2 



4DJ 



(30) 




Figure 7: Diagramatic representation of incoherent scattering 
processes resulting in induced losses and light shifts. 



where Di = cjo ^ and £'2=^0" ^2 represent the de- 
tunings of lasers uji and uj2 from the resonance frequency 
uoq of the atomic transition between the dissociation Hm- 
its of the ground and excited potentials. In addition. 



/3I = 



^'-^2v',UJ2 



^'^2v',ui 



Iv' ,UJ\ Iv' 



Iv' ,CJi 



2V2 



h,v'l2,v' 



(31) 



(32) 



(33) 



(34) 



Dj + i-fA/2 , 



(35) 



where ja = 7/2 is the atomic spontaneous decay rate, 

(A) (A) 

and fl\ and are the atomic Rabi frequencies which 



we take n^"^^ = 

We now introduce real frequency shift and loss coeffi- 
cients, so that = I3i -\- ia/2 , jS^ = Pi + 1^2/'^ and 
tj^^ = IJxi+iTi/2 . Note that, in general, x* 7^ x' ■ This 
means that the parametric terms are not population- 
preserving, and can provide a STIRAP-type of loss re- 
duction even for the non-resonant vibrational levels, pro- 
vided the coupHng is STIRAP-Hke. 

The coefficients a, F^, (3i, Uu and x are obtained by ex- 
plicitly treating all other levels in the excited potential, 
adjacent to |3), followed by the procedure of adiabatic 
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elimination. The coefficients also include the contribu- 
tions from Raman type of couplings |1) \v') by the uji- 
laser and |2) ^ \v') by a;2-laser. In principle, these addi- 
tional Raman couplings could be treated exactly like the 
primary STIRAP transition via |3), i.e. taking place via 
the dark-state interference effect, except that the transi- 
tions have much larger one-photon detuning. This would 
require an adiabaticity condition of the form of Eq. ( [T^ ) 
that includes the one-photon detuning, implying that a 
larger value of the product n^f-^-^'^''T is needed. However, 
^(e//,o)j, made arbitrarily large, as we dis- 

cussed earher. Therefore our approach is to treat these 
extra Raman couplings as loss and dephasing processes, 
rather than to include them into the adiabatic passage 
scheme. The contribution of these coupHngs to the ef- 
fective atom-molecule conversion rate, described by x, is 
negligibly small compared to the conversion rate due to 
the parimary Raman transition via |3). 

The relevant transitions that stand behind these co- 
efficients are illustrated in Fig. |^. For example, the 
coefficient a describes the process of atomic absorption 
from either of the two Raman lasers, that incoherently 
produce excited atoms followed by spontaneous emission 
loss. The coefficient Fi is due to ordinary photoassocia- 
tion when pairs of atoms from the condensate are trans- 
ferred (again by either of the two Raman lasers) into 
an excited molecular state which can then spontaneously 
dissociate into a pair of hot (non-condensed) atoms. The 
effective rate of this non-linear loss is Firii. Finally, the 
coefficient T2 describes the loss of molecules due to spon- 
taneous Raman scattering of laser photons. This pro- 
duces molecules predominantly in ro-vibrational levels 
other than the one targetted by the stimulated Raman 
transition. This term may also describe some scattering 
which is elastic in the sense that the molecules return to 
the targeted state, but with an increased kinetic energy 
due to photon recoil, that will still remove the molecules 
from the condensate. 

The summations in the expressions for F^, (32, Uu and 
X are carried out over all the excited levels v' except the 
resonant level |3) which itself participates in STIRAP, 
rather than being adiabatically ehminated. The effects 
of losses and light shifts due to the cross-coupHngs to 
the level |3) itself are implicitly described by the last 
terms in the rhs of Eqs. (p7|)-(p9|). Subsequently, we will 
estimate the combined effects of all levels in which case 
the contribution of the level |3) is estimated by similar 
terms to the ones included in Ti, P2, Uu, except that the 
detunings Aii,/,^^ and A2v/,uji are replaced by Ai^^^^^ ^-^d 
^23,1^1, respectively, and li^v' are replaced by li^^. 

Separating out the time dependences of the two Rabi 
frequencies, the above coefficients can be rewritten as: 



(36) 



= ^(i)e-2(*-t,)VT= +^(2)e-2(t-*.)VT= 
X = Xoe~ 

where the peak values are, respectively: 



4 



^1 - J 2^,, 



{A,0} 



D, 



7 y^' 



(el.O) 



A21,', 



(^,0)|2 



4A 



P'2'- 



Xo 



Xo 



^(e',0)^(ei,0)* 



Iv ,a;i 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



p(l)g-2(t-ti)VT= 



■2it-t2f/T^ 



The reason for this separation is that the two terms in 
each coefficient act during different time intervals, corre- 
sponding to the first and the second pulse in STIRAP. 
Accordingly, one has to distinguish their disruptive effect 

during the duration of the corresponding pulses. For ex- 

(2) 

ample the molecule loss term Fj acts during the first 
Raman pulse (of frequency 0^2) when the molecular field 

(2) 

is not populated yet. As a result, the coefficient Fj is 
not so disruptive. On the other hand, the molecule loss 
term Fj^^ is much more important since it acts during 
the second Raman pulse (with frequency uji) when the 
population of the molecular condensate becomes high. If 
the value of F2^^ is too large, one can easily lose all this 



(37) population during the wi-pulse. 



11 



In order that the radiative losses and dephasing due 
to light shifts be negligible over the duration of STIRAP, 
the time scales associated with the coefficients a, Fini, 
r2 , and the induced relative phases must be much larger 
than the duration of pulses in STIRAP, i.e. 



a 



(48) 



(49) 



.(2) , 7 



"2 ^1,3 



ni 



.(1) ^ 7 



"l -'2,3 



> T, (50) 



(51) 



E 2 




160 



0.3 



0.2 



0.1 



170 180 

v' 



190 



(b) 




50 



100 
v' 



150 



200 



.(2) 



(52) Figure 8: Free-bound (a) and bound-bound (b) Franck- 
Condon overlap integrals, Ji,,,' and /„^35,j,/, as a function of 
the vibrational quamtum number v' . 



Pi'' 



4A 



23,cji 



2(3\ 



(1) 



(53) 



Pi'' 



(2) 



(54) 



are minimized. This is done using the results of calcula- 
tion of the dipole matrix elements, energy eigenval- 
ues, and the Franck-Condon overlap integrals for model 
potentials that closely approximate the ^^Rb2 ground 
^ J2u potential and the 0~ symmetry excited potential. 
The calculation treats 205 ro- vibrational levels in the ex- 
cited potential (which we label by v' = 0, 1, 2, 204) 
and 39 levels {v = 0, 1, 2, 38) in the ground potential. 

Within such a large range of target levels that the Ra- 
man transitions can be tuned to, several possibilities can 
be readily eliminated to simplify the search. For exam- 
ple, Raman transitions via one of the highly excited levels 
{v' > 190) will suffer from large values of the atomic loss 
coefficient a since the detunings Di will be small, and 
as a result the condition [a^*-'] 3> T will not be satis- 
fied. On the other hand, transitions via low excited states 
{v' < 160) will have very small values of the free-bound 
Franck-Condon overlap integral, Ii^v' ^ 0.5 x lO^^'' m^^^. 
This in turn will result in small effective Rabi frequency 
il'f-^^'^^ (using reasonable values of the intensity of the 
\ and f3y. Therefore these additional contributions laser uji and the density ni), so that the adiabaticity 
must be included in the above conditions, as they play an ,o„dition T » 1 is not satisfied, 

important role for correct estimates of the overall degree 

of disruption due to incoherent coupHngs. In general, the behavior of the coefficients Ty' , F^'^ 

Our goal now is in performing a realistic analysis of the ""-^ '^^^^ 
above coefficients for the ^''Rb2 molecule under consider- 
ation, and in finding appropriate target levels for the Ra- 
man transitions in STIRAP so that the disruptive effects 



This will guarantee that the pulses are switched off before 
the losses and dephasings can have their disruptive effect. 
The infiuence of the nonlinear phase shift due to u[i can 

be ignored simply on the grounds of t/j;/ ^ t/ii which is 
the case we encounter in our analysis. 

(2) 

In the above conditions involving the coefficients T\ , 

Fj^-*, and l3i^\ we have included additional terms which 
are the contributions from the incoherent cross-couplings 
1 <-*■ 3 by the laser 102 and 2 ^ 3 by the laser uji. As we 
mentioned earlier, these processes are treated explicitly 
by the last terms in the rhs of Eqs. (p7|)-(p9|). However, 

their overall effect can be described by expressions sim- 

(2) 

ilar to the corresponding terms in the coefficients T\ , 



,(1) 



and /Sj is not of a trivial character. Each particular 
choice of the final state \v) in the ground potential which 
we designate as |2) will result in different sets of the de- 
tunings of the Raman lasers with respect to couplings 
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51.43 s"^ 




0.5398 


r(i) 


3.010 X 10"^'' mVs 


r(2) 
^ 1 


3.014 X 10-^^ mVs 


p(i) 

^ 2 


466.7 


r(2) 

^ 2 


4.645 s"^ 




2.948 X 10® s"^ 




3.020 X 10 s 




5.652 X 10'^ s"^ 




5.969 X 10'' s"^ 


(7/4) 




'/|Ai3,.,r 


7.96 X 10~^^ m^/s 


(7/4) 




'/|A23,^J' 


375.9 


|n(^'-°)j2,3|V(4A23,.J 


2.57 X 10^ s"^ 


x'o ~ Xo 


-9.25 X IQ-'' m^/^/s 



(a) 



Table IV: Calculated values of the loss and light shift coeffi- 
cients for Raman transitions tuned to the ground t; = 35 and 
excited u' = 177 levels. 



to different excited levels. Further complications emerge 
from the oscillatory behavior of the Franck-Condon over- 
lap integrals, as shown in Fig. ||. This means that 
the contribution of the levels nearest to |3), having the 
smallest detunings, may not necessarily give the leading 
term in the sums over v' , since the further detuned levels 
may have larger Franck-Condon overlaps thus resulting 
in comparable contributions to the coefficients. In addi- 
tion, different terms in the expression for will depend 
on the sign of the respective detuning, so that they may 

(i) 

add up into either a positive or negative value of ■ 

Thus our analysis consists of a calculation of all the 
above coefficients for different levels \v) and \v'), and 
subsequent identification of a n op timum target level that 
satisfies the conditions (^-(^) as closely as possible. 
For each level \v) = |2) in the ground potential we scan 
the Raman transitions through different levels \v') in the 
excited potential, treating this as |3) and carrying out 
the summations over the remaining levels. 

The calculation is done for Sl^'^''*''' = 10^" s~^ and 
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Figure 9: (a) Optimum effective Rabi frequencies for STI- 
RAP in the multi-level model, (b) Resulting densities in the 
uniform condensates case. 



tential bound by -8.05 GHz or -5.058 x lO^o s"!. The 
binding energy of the level u = 35 sets up the frequency 
difference a;2 — wi = 5.058 x 10^° s~^, and the values of all 
relevant detunings. The nearby levels around w' = 177 in 
the excited potential are separated by about 4.2 x 10^^ 
s^^, which is much larger than W2 — wi. For this arrange- 
ment of the target levels, the resonant Franck-Condon 
overlap integrals are equal to /i_3 = 1.05 x lO"^"* vo?^^ 
and /2,3 = 0.0228, so that the effective peak Rabi fre- 



,(e//,0) 



= 2.18 X 10'' 



(ei,0) 



2 — 10^ s~ . To simplify the analysis we consider 
the case of Ai =0 and (5 = 0, i.e. we only treat cases 
when the primary Raman transition is resonant. Once 
the optimum level is identified, we can subsequently fine 
tune the two-photon detuning 5 for optimum conversion. 
Provided that the order of magnitude of the optimum 5 is 
\5\ ~ 10** s~^ as found in the examples of sub-sections IV. 
B and C, this subsequent fine tuning will have no effect on 
the results of calculation of the loss and light shift coeffi- 
cients, since these involve detunings i^iv' ,u:j that typically 
have much larger magnitudes, |Ai^/,i^J ^ \5\. 

The most favorable case that we find corresponds to 
tuning the Raman transitions to v' = 177 in the excited 
potential bound by -22.23 cm'^ or -4.1903 x lO^^ g-^i 
from the dissociation limit, and w = 35 in the ground po- large and the adiabaticity condition is maintained. 



for 

-1 



quencies are equal to: VL\ 

ni(0) = 4.3 x 10^0 m-3, and d^^^'^'^ = 2.28 x 10^ s 
The resulting values of calculated loss and light shift co- 
efficients are given in Table IV . 

An important factor in minimizing the most significant 
loss coefficient Fj^^ and its "cousin" term in Eq. (^ ) 
is the relatively small value of /2,3 and of the Franck- 
Condon overlaps of the closest nearby levels. While 
this is favorable for the undesired loss terms, the small 
value of /2,3 also affects the strength of the bound-bound 

coupling of the primary Raman transition, i}^-^^'^^ — 
j-^(ei,o) j^^^ which must be kept large. However, the small 
-f2,3-value is compensated here by a strong bare electronic 
Rabi frequency fl'^^'^^ = 10^ 



1 so that ni""^^'"^ is still 
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Using the above parameter values and simulating the 
STIRAP equations with the additional terms given by 
Eqs. (p7|)-(p9|), results in about 42% conversion efficiency 
{r] — 0.42), which is rather high and is an encouraging 
result. In this simulation, the estimated optimum pulse 
duration T, the pulse delay coefficient a, and the two- 
photon detuning S, were taken as follows: T = 1.7 xlO~^ 



(a) 



1.53, and (5 = 4 X 10* s" 



The effective Rabi 



frequencies and the particle number densities for this cal- 
culation are given in Figs. || (a) and (b). We note that 
these results are essentially identical whether the exact 
complex form or the far-off-resonant approximations to 
the additional coefficients are employed. In other words, 
any coherent cancellation process via STIRAP in the ex- 
tra vibrational levels, is neghgible compared to the inco- 
herent losses via non-STIRAP excitation of these levels. 



B. Non-uniform condensates 

The final step in our analysis is to include the trap po- 
tential and the kinetic energy terms and simulate the full 
set of coupled inhomogeneous mean-field equations (^) in 
three space dimensions. The initial state in these simu- 
lations is a pure atomic BEC, with no molecules present, 
as given by the standard steady-state Gross-Pitaevskii 
equation in a trap. We consider spherically symmetric 
trap potentials X^(x) = Ei + (mi/2)a;||xp and choose 
the trap oscillation frequencies uii equal to each other: 
uJi/2Tr = 100 Hz (i = 1,2,3). Including these terms, we 
simulate Eqs. (g) assuming that the initial peak den- 
sity of the atomic BEC is ni(x = 0,t = 0) = 4.3 x 10^° 



This corresponds to the total 
iVi = / d3x|^/>(x,0)p equal to 



m at the trap center, 
initial number of atoms 
Ni = 5x 10^ 

We note that the characteristic time scale associated 
with the trap potential can be estimated as = l/w ~ 
1.6 X 10~^ s. The time scale associated with the ki- 
netic en ergy term is estimated from the 'healing' length 
Ih ^ y/htfi/mi which corresponds to a healing time 
scale of th ■ Under adiabatic conversion of an equilibrium 
(Thomas-Fermi like) atomic BEC, this healing time scale 
coincides with the dephasing time tph associated with the 
mean field energy potential and discussed before. Both 
these time scales (i^j and th) are longer than the dura- 
tion of pulses in STIRAP we employed earher. Therefore, 
addition of these terms can not dramatically change the 
results and conclusions obtained above for uniform con- 
densates. 

To show this we first consider the idealized three-level 
model with parameter values given in Tables | and ||, 
i.e. the case of equal and relatively strong effective Rabi 
frequencies. The results of simulations of Eqs. (||) are 
shown in Figs. |l^ and 11, where we see a rather high 

0.91. This should be compared 



|10| and 

conversion efficiency 77 ~ 

with the uniform condensates result of Fig. |[ 

Next, we consider the case of more reaHstic parameter 
values (lower Rabi frequencies) , and include the effects of 




4 0-6 
n 0.2 "-^ 
t (ms) 



Figure 10: Densities ni(x,t) = \ipi(x,t)f of the atomic (a) 
and molecular (b) condensates in a trap as a function of time 
t and the radial distance r = lx| from the trap center. The 
applied pulses are as in Fig. g, and other parameter values 
are given in Tables Q and |l| The result is for the idealized 
three-level model. 



incoherent radiative losses and dephasings as described in 
the previous sub-section. Using the calculated values of 
all relevant parameters for the Raman transitions tuned 
to v' = 177 and v = 35, and adding the trap potential 
and kinetic energy terms to the earlier Eqs. ([2^)-(29), we 



simulate the process of STIRAP governed by the result- 
ing full set of mean-field equations in three dimensions: 



dt 



■AS"-*, 

a 



V2 



ri^3, (55) 



dip2 (x,t) 
dt 



-V'3 



r2 



,11^12* 



^3, (56) 
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Figure 11: Integrated occupation numbers Ni{t) = 
J dyi\^i{x.,t)f of the atomic (solid line) and the molecular 
(dashed line) fields as a function of time t, for the parameter 
values of Fig. 
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Figure 12: Same as in Fig. |ll]but for the parameter values of 
the complete multi-level model and a three-dimetional trap 
geometry. The resulting conversion efficiency can be com- 
pared with the uniform condensates case of Fig. ^ (b). 



dt 



2^/2 

^{"^^h 3 ■ 



(57) 



We switch on the sequence of two Raman pulses, as 
shown in Fig. 9, corresponding to: = 2.18 x 10® 



T = 1.7 xlO- 



s~\ n^^''^^'°'> = 2.28 X 10^ s % J = i.Y xiu ' s 
and a = 1.53. The values of T and a, together with 
the choice of the two-photon detuning S — A x 10'' s~^, 
correspond to the optimized set of parameters as in Table 
III and are obtained for v' — 177 and v = 35 levels 



using il^ 



(el,Q) 



10 



10 



and Q 



(ei,0) 



10^ 



The 



corresponding values of all other relevant coefficients are 
given in Table IV and while the spontaneous decay 
rate is 7 = 7.4 x 10^ s^^ , as before. The optimum value of 
5 accounts for nonzero values of the light shift coefficients 
Pi , so that the actual detuning to be optimized using the 
Table [II is the overall effective detuning S = S+f32 — 2/3i . 
The results of simulations with this set of parameter 



values are given in Fig. |l^ where we plot the total oc- 
cupation numbers in the atomic and molecular BECs as 
a function of time. The obtained conversion efficiency is 
77 ~ 0.32. This is 10% lower than the efficiency in the 
corresponding homogeneous case of Fig. but still is a 
rather encouraging result given the fact that 32% conver- 
sion of about 5 X 10^ atoms would give a molecular BEC 
with the total of 8 x 10'' molecules, with a peak density 
of about lO^o m-3. 



VI. SUMMARY 

To summarize, STIRAP is a potential route towards 
the coherent conversion of an atomic to a molecular BEC. 
This process involves stimulated emission of molecules, 
and is very different to normal chemical kinetics. As 
such, it is a form of 'superchemistry' STIRAP in 
an atomic BEC can be treated in a very similar way to 
normal STIRAP, by introducing an effective Rabi fre- 
quency for the first photo-association transition. A po- 
tential difficulty is the relatively low values of the effective 
Rabi frequency in the first or photo-association transi- 
tion. This can be regarded as physically due to the low 
densities of atoms in a typical weakly interacting BEC, 
compared with atoms in a molecule. This means that 
the corresponding Franck-Condon coefficient, after mul- 
tiplying by the relevant BEC amplitude, has a very small 
value. Thus, the laser intensities required may be quite 
high, in order to obtain Rabi frequencies comparable to 
those used in atomic transitions. 

This by itself is not critical, since deep in the adiabatic 
limit it is normally permissible to use low Rabi frequen- 
cies, as long as the associated time-scales are long enough. 
From this point of view, the use of STIRAP, and conse- 
quent reduction of spontaneous emission, is a physically 
sensible idea. However, including s-wave scattering or 
mean-field processes into the model for STIRAP sets up 
a characteristic two-photon dephasing time scale, so that 
the pulse durations in STIRAP have to be shorter than 
a certain critical value. Short pulse durations necessar- 
ily involve high values for the effective Rabi frequencies 



in the usual symmetric case of fli 



(e//,0) 



(e//,0) 



, thus 



requiring a very high laser power for the first transition 
- if the adiabaticity condition is to be maintained. 

In order to ease this demanding requirement and be 
able to achieve highest possible conversion efficiency at 
smaller total laser power, we propose to use an off- 
resonance operation (thus cancelling part of the mean- 
field detuning effect) and effective Rabi frequencies 



(e//,0) 



and n 



(e//,0) 



of different magnitudes. This has 



the effect of dramatically increasing coherent molecule 
production, in a physically accessible regime of mod- 
erate laser intensity. Further improvements may be 
possible by tailoring the input pulse frequencies to the 
time-dependent two-photon detuning, caused by inter- 
atomic and intermolecular scattering. We have carried 
out mean-field calculations in three dimenisons to verify 
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that trap inhomogeneity should not have adverse eflects 
on the STIRAP process. 

Finally, we stress the importance of radiative losses 
and dephasing due to incoherent processes that occur 
during STIRAP. These processes are usually assumed 
to be negligible in ordinary STIRAP between purely 
atomic or molecular states. In the present case of coupled 
atomic/molecular BECs, this assumption can not be eas- 
ily justified since the free-bound transition typically in- 
volves a relatively low effective Rabi frequency, which ne- 
cessitates long pulse durations. Instead, we find that the 
incoherent couplings can be rather destructive unless spe- 
cial care is taken to minimize their effect. This involves 
detailed knowledge of the structure of the free-bound and 



bound-bound transitions and subsequent identification of 
optimum target levels in STIRAP, so that the overall con- 
version efficiency remains comparable to the predictions 
of the simplified three-level model. 
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